Nano-Engineering Defect Structures on Graphene 
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We present a new way of nano-engineering graphene using defect domains. These regions have 
ring structures that depart from the usual honeycomb lattice, though each carbon atom still has 
three nearest neighbors. A set of stable domain structures is identified using density functional 
theory (DFT), including blisters, ridges, ribbons, and metacrystals. All such structures are made 
solely out of carbon; the smallest encompasses just 16 atoms. Blisters, ridges and metacrystals rise 
up out of the sheet, while ribbons remain flat. In the vicinity of vacancies, the reaction barriers to 
formation are sufficiently low that such defects could be synthesized through the thermally activated 
restructuring of coalesced adatoms. 



Carbon is a fundamental material for nano- 
engineering. Three dimensional graphite, one- 
dimensional nanotubes, zero-dimensional buckyballs, 
and now two-dimensional graphene are all being investi- 
gated intensely to this end. The latter, recently created 
experimentally for the first time [H, 0] , is the first stable 
two-dimensional solid state lattice material. Fundamen- 
tal aspects of graphene include the room-temperature 
quantum Hall effect and an effective description by a 
Dirac equation |3|. Graphene has significant applications 
in electronics [JlU Q ! and is even predicted to replace 
silicon in future solid state devices Many of these 
applications require cutting sheets of graphene precisely 
in order to localize charge, among other effects. We 
propose an alternative method of patterning graphene 
without cuts of any kind, which can be adapted to 
many applications. The key to this method is controlled 
placement of groups of defects, called defect domains. 

In this Letter, we demonstrate that a complex, stable 
landscape of defect domains can be monolithically nano- 
engineered from graphene to make arbitrary structures. 
An example of a defect domain which stands up out of 
the graphene plane is shown in Fig.[lja)-(b). Such defect 
blisters in graphene can be generated through the coales- 
cence of adatoms in the vicinity of vacancies (Fig. [He)). 
Defect ridges can be constructed as an alignment of blis- 
ters and their two-dimensional patterning results in de- 
fect metacrystals. Planar, defect ribbons may also be 
synthesized as shown in Fig. QJd). These structures, in 
particular, may be viewed as being spliccd-in pieces of 
a carbon allotrope known as haeckelite Q • Figure Ee)- 
(f) show how haeckelite symmetry considerations may 
be used to predict energetically favorable defect domain 
geometries. We use density functional theory (DFT) to 
study stability, scalability, and other key physical proper- 
ties of these materials. DFT has been applied successfully 
to magnetism in graphene nano-islands [H, the effect of 
substrates Q , and transport in doped graphene nanorib 
bons [10| ; DFT is also used for carbon nanotubes 111, [ 
and other carbon structures [HI]. 

Defects have been observed in graphene [H, [3| and 
are expected to play a key role in the functional prop- 
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FIG. 1: Defect Domains in Graphene. (a)-(b) Top and side 
view of a single defect blister of 24 carbon atoms with a height 
of 1.9 A. (c) A Stone- Wales defect used to synthesize defect 
domains, (d) A defect domain ribbon embedded in graphene. 
(e)-(f) Haeckelite, a carbon allotrope used to synthesize defect 
domains. In yellow is shown a graphene overlay for visualiza- 
tion. 



erties of working materials fl5j | . Among the simplest 
and more easily fabricated of these is the Stone- Wales 
defect, wherein the rotation of a single pair of carbon 
atoms creates adj acent pairs of pentagonal and hep- 
tagonal rings |16j, as shown in Fig. [ljc). Such de- 
fects can be introduced intentionally using electron radia- 
tion and imaged using Transmission Electron Microscopy 
(TEM) Q. Electron radiation also causes the forma- 
tion of adatom/ vacancy pairs which subsequently sep- 
arate and move across the graphene lattice DFT 
analysis ljj estimates the adatom migration barrier to 
be 0.45 eV with a jump frequency of (3.7±0.7) x 10 12 s~\ 
while the barrier to vacancy migration is estimated to be 
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1.7 eV with a jump frequency on the same order as for 
adatoms 18j. Electron radiation can therefore be used 



to generate adatoms and vacancies which have been ob- 
served to rearrange themselves to form more complex de- 
fect structures [14j . 

Domains over which similar defect structures are peri- 
odically replicated are known as haeckelite p}, as shown 
in Fig. WIS)- Tight binding molecular dynamics stud- 
ies suggest that small regions of these carbon allotropes 
may be induced within graphene through the coales- 
cence of four single vacancies [19j . In principle, entire 
sheets and tubes of haeckelite can be synthesized which, 
even without experimental realization, have already gen- 
erated potential applications [23]. We focus on the sim- 
pler notion of localized defect structures which are mini- 
mally extended beyond Stone- Wales defects. We initially 
treat defects composed of ribbons and patches of haecke- 
lite and then consider a new fundamental building block 
which we term an Inverse Stone- Wales defect. 

All calculations were performed with the real-space, 
numerical atomic orbital, DFT code, DMOL 2l|. A 
norm conserving, spin unrestricted, semi-core pseudopo- 
tential approach was employed with electron exchange 
and correlation accounted for using the Perdew-Wang 
generalized gradient approximation (GGA) [2^]. Periodic 
boundary conditions were employed and vacuum slabs 
were used to isolate the replicated graphene layers. As a 
check on the method, the ground state energy of Cgo, i.e., 
a buckyball, was estimated to be 384 meV/atom above 
that of graphene, consistent with a literature value of 
380 meV/atom [231 ] . Likewise, a single Stone- Wales de- 
fect was estimated to have a formation energy of 5.08 eV 
when embedded within a 144-atom graphene supercell; 
this compares well with an estimate of 4.8 eV from the 
literature [2lj ]. 

We begin with defect structures constructed from 
patches of haeckelite Q. The Hg g,7 variant with hexag- 
onal symmetry was chosen over rectangular and oblique 
allotropes. This is the most stable of the three variants. 
We find a ground state energy estimate for H 5:6i 7 of 229 
meV/atom above graphene. This is close to the liter- 
ature value of 246 meV/atom and notably lower than 
that of Ceo [Hj] . H 5 , 6j 7 is 3% less dense than graphene: 
it has 0.369 atoms/A 2 , compared to 0.380 atoms/A 2 for 
graphene. Haeckelite can therefore be formed by the 
removal of carbon atoms or in settings in which the 
graphene sheet can expand. For instance, a 20-atom pe- 
riodic cell of this carbon allotrope can be formed through 
the addition of two carbon atoms to a Stone- Wales defect 
followed by lateral dilation. The energy barrier is esti- 
mated to be 1.8 eV per cell. Such restructuring should 
be aided by straining the graphene. DFT calculations 
support this and indicate that the stability exchange oc- 
curs at a linear strain of ~0.07, as shown in Fig. O The 
analysis also determined uni-axial elastic constants for 
graphene and haeckelite of 1.13 GPa and 1.20 GPa, re- 
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FIG. 2: Mechanical response (pure dilation) of graphene, 
haeckelite, and a metacrystal of haeckelite-based defect do- 
main blisters. The points represent actual DFT data; the 
curves are a guide to the eye. 



spectively. The first number is consistent with experi- 
mental measurement of 1.06 GPa 25[, while the second 
indicates that haeckelite is 6% stiffer than graphene. 

Two types of approximately coherent interfaces be- 
tween graphene and H5 ^,7 can be identified by inspection 
and are shown in Fig. [ije)-(f). A linear haeckelite strain 
of 4.7% gives averaged coherence along the symmetry 
boundaries shown in Fig. [He), while triangular patches 
of haeckelite are approximately coherent without apply- 
ing any dilation, as illustrated in Fig. WIS)- With these 
lines as a guide, haeckelite can be spliced into a graphene 
sheet and vice-versa. A strip of haeckelite spliced into 
graphene along adjacent, linear coherency lines has a pla- 
nar ground state, as shown in Fig.QTd), where the dashed 
lines indicate the periodicity of our simulation. On the 
other hand, splicing a triangular tile from Fig. [Iff) into 
graphene causes a ground state blister to form, as shown 
in Fig. Dla)-(b). The converse, a triangular graphene 
patch within haeckelite, results in a planar structure. 

The concept of a defect blister creates a natural way 
to extend graphene out of the plane. A formation en- 
ergy of 6.07 eV was determined for the blister illustrated 
in Fig. QTa)-(b), using a 202-atom supercell. This is ap- 
proximately 1 eV higher than the Stone- Wales defect. 
A Hessian analysis [26| on a smaller supercell indicates 
that such haeckelite-based defect blisters are linearly sta- 
ble. This is consistent with the stability observed in a 
room temperature quantum molecular dynamics simula- 
tion, which we also performed. 

We now turn to an even simpler defect domain struc- 
ture not based on Haeckelite, and explicitly describe its 
reaction pathway. While haeckelite-based defect blisters 
amount to the substitution of a pair of atoms for a hexag- 
onal ring (Fig. [Hf)), a narrower structure can be formed 
by inserting two atoms on opposing faces of an exist- 
ing hexagon as shown in Fig. OHa)-(b). This results in 
pairs of pentagons and heptagons that we refer to as an 
Inverse Stone- Wales defect (compare to Fig.QJc)). A 16- 




FIG. 3: (a)-(b) A single Inverse Stone-Wales defect, (c)-(f) 
Multiple defects can be aligned to form ridged contours which 
can be straight, curved, or even closed. 



atom blister results with a footprint of 12.2 x 7.4 Aand 
a height of 2.1 A. A formation energy of 6.20 eV was 
determined using a 200-atom supercell, nearly the same 
energy as that for the larger haeckelite blister. A series 
of Inverse Stone- Wales blisters can be aligned to form a 
corrugated ridge of nearly arbitrary contour (Fig. G3c)- 
(d)), a monolithic analog of carbon nanotubes formed 
from Ceo molecules 27]. Such structures may hold util- 
ity in guiding charge transport on graphene. Further, the 
creation of closed contours results in extended blisters of 
pure graphene enclosed by a ring of Inverse Stone- Wales 
defects. (Fig. |U(e)-(f)). 

The synthesis of defect structures may be facilitated by 
electron radiation and the resulting collection of adatoms 
and vacancies [14] - The formation energy for a single 
vacancy was calculated to__be 7.63 eV, consistent with 
experimental estimates 



28], and therefore has the po- 



tential to provide a large driving force for restructuring 
events within its horizon. The formation energy of a va- 
cancy/adatom pair, the planar analog of a Frenkel defect, 
is 14.13 eV. Di-vacancy formation energy was calculated 
to be 8.08 eV. In addition, re-structuring within such a 
low dimensional system is more easily carried out than 
in bulk since vacancies facilitate significant low energy 
distortions of the local lattice. To explore this, a hybrid 
linear synchronous transit/quadratic synchronous tran- 
sit (LST/QST) transition state search algorithm [29[ was 
used to construct a reaction pathway for the approach of 
a single adatom to a di-vacancy (Fig. |4j. Adatoms hop 
between adjacent bridge sites, and a barrier of 0.52 eV 
was identified for jumps in pure graphene, consistent with 
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FIG. 4: (a)-(b) Reaction path for a single adatom moving to 
the periphery of a di-vacancy. (c) Experimentally observed 
adatom positions [31 . (d)-(e) Creation of an Inverse Stone- 
Wales defect on the periphery of a di-vacancy by one adatom 
hopping within the vicinity of another, (f) Final structure 
with blister shown in dark gray and di-vacancy identified 
within a dashed ellipse. Energies in both (a) and (f) are 
relative to the lowest point plotted. 



an estimate of 0.45 eV obtained elsewhere [Uj. How- 
ever, the di-vacancy attracts adatoms to positions that 
are one hop removed from its periphery (position 3 in 
Fig. Htb)), and the associated barriers can be overcome 
by thermal fluctuations at modest temperatures. These 
are precisely the bridge sites at which adatoms are ob- 
served experimentally (Fig. Etc)) 14j. The reaction en- 
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ergies in Fig. (Ha) are E\, 
eV, and E^ 4 = 2.20 eV, while the reaction barriers are 
E\ 2 = 0.14 eV, B| 3 = 1.55 eV, and E b 3i = 2.49 eV. 

Two adatoms are required to synthesize an Inverse 
Stone- Wales blister. So a second transition state analy- 
sis was carried out wherein one adatom jumped towards 
another located at a low energy site on the periphery 
of the di-vacancy, as shown in Fig. IUe)-(f). As indi- 
cated in the figure, the reaction barriers can be overcome 
by thermal fluctuations at modest temperatures and re- 
sult in a blister. The reaction energies are E\ 2 = —0.40 
eV and K 
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-7.10 eV, while the reaction barriers are 
E\ 2 = 0.39 eV and E^ — 0.34 eV. Subsequent vacancy 
migration or elimination would leave the stand-alone blis- 
ter shown in Fig. [3fa)-(b). 

Localized defect domains are of interest in their own 
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FIG. 5: (a)-(b) Top and side view of a metacrystal formed 
from haeckelite-based defect-domain blisters, (c) A metacrys- 
tal defect: a less closely packed metacrystal contains an 
"atom" of negative polarity, a blister pushed down rather than 
up. 



right, but patterned arrangements of them may en- 
dow graphene sheets with novel metacrystalline proper- 
ties. Within such a paradigm, blisters represent meta- 
atoms with either up or down polarity. A mono-sized, 
mono-polarity, ground state, metacrystal is shown in 
Fig. E^a)-(b). The associated mechanical response given 
in Fig. [2] indicates that it has a density greater than 
that of graphene. The ground state energy is only 274 
meV/atom above that of graphene and is surprisingly 
close to the ground state energy of pure haeckelite, which 
is 229 meV/atom above graphene. This suggests that the 
metacrystals are more energetically stable than Ceo; the 
latter has an analogous energy of 384 meV/atom. The 
elastic stiffness of the metacrystal is 1.08 GPa in uni-axial 
tension, somewhat less that either graphene or haeckelite. 
One possibility for the experimental realization of defect 
metacrystals in graphene is interference patterns of elec- 
trons in the near field of a diffraction grating [3(|. Such 
a system has been demonstrated at the 100 nm scale al- 
ready and has the potential to be reduced to the 10 nm 
scale. Metacrystals can also be made with larger blisters 
or even blisters of alternating size. Thus one can consider 
a "periodic table" of meta-atoms, creating a variety of 2D 
metacrystals of arbitrary lattice and crystal structure. 

In conclusion, we have described a new method to 
achieve monolithic nano-engineering of graphene via de- 
fect domains. These structures can take the form of blis- 
ters, ribbons, and ridges. We showed that such objects 
are linearly stable and can be arranged in arbitrary pat- 
terns, leading to a metacrystal. The smallest defect do- 
main, referred to as an Inverse Stone- Wales defect, con- 
sists of only 16 atoms, and the associated energy is ap- 
proximately 1 eV higher than a simple Stone- Wales de- 



fect. Di- vacancies were shown to attract adatoms to their 
periphery and a thermally activated path was identified 
for blister synthesis. Defect domains may offer techno- 
logical applications associated with the confinement and 
transport of charge, as we will investigate in future work. 

We acknowledge useful discussions with James 
Bernard, David Wood, and David Wu. LDC was sup- 
ported by the National Science Foundation under Grant 
PHY-0547845 as part of the NSF Career program. 
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